Phase separation in binary colloids with charge asymmetry 



<N 

o 

<N 



m 



o 



O 

o 



> 

CO 

p 

en 

o 

(N 



't 



Koki Yoshizawa^, Nao Wakabayashi^ , Masakatsu Yonese^, Junpei Yamanaka^ and C. Patrick Royall^ 

^Faculty of Pharmaceutical Sciences, Nagoya City University, 

3-1 Tanahe, Mizuho, Nagoya, Aichi 467-8603, Japan, and 

^School of Chemistry, University of Bristol, Bristol, BS8 ITS, UK}^ 

(Dated: February 24, 2013) 

We report that binary dispersions of like-charged colloidal particles with large charge asymmetry 
but similar size exhibit phase separation into crystal and fluid phases under very low salt conditions. 
This is unexpected because the effective colloid-colloid pair interactions are accurately described by 
a Yukawa model which is stable to demixing. We show that colloid-ion interactions provide an 
energetic driving force for phase separation, which is initiated by crystallization of one species. 
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I. INTRODUCTION 

We encounter suspensions of charged colloids every 
day: their applications range from cosmetics to coatings. 
Since they are well described by a soft repulsive inter- 
action, charged colloids also form an important model 
system, and behave in a similar way to materials as di- 
verse as dusty plasmas [U [2] and metals [3^. Further- 
more, colloidal particles in aqueous suspension form a 
basic model of biological systems such as proteins ^ . In 
all these examples, many-body effects (deviations from 
density independent pairwise interactions) can be impor- 
tant: many-body effects underly the phase separation we 
observe here. 

Charged colloids in suspension interact with small ions 
in the form of added salt, and counterions which balance 
the charge on the colloids. This system can in princi- 
ple be described with Coulomb interactions between the 
small ions and colloids, with an additional hard core term 
to account for the size of the latter. However, usually the 
number of interacting species is too large to deal with. 
The degrees of freedom of the small ions may be inte- 
grated out, which leads to an effective colloid Hamil- 

tonian W = Wi{ps,pc;pk) + W2{ps, Pc] pkA^}) + ••• 
where the one-body terms Wi do not depend on the col- 
loid coordinates {R}. Here the counterion number den- 
sity pc = Pion — "^Ps where pion is the number density of 
(monovalent) ions, ps is the number density of salt ion 
pairs and pk is the number density of colloid species k [5 . 
This effective Hamiltonian then depends on state point 
and is not in general pairwise additive. 

Often, one considers the effective colloid-colloid terms 
W2 only, neglecting higher-order and one-body terms, as 
the latter are independent of colloid coordinates. The 
two-body terms involve an effective interaction between 
the colloids which can be described by the Derjaguin, 
Landau, Verwey and Overbeek (DLVO) model [6 . This 
assumes that the electrostatic potential is sufficiently 
weak that linear Poisson-Boltzmann theory can be ap- 



plied. Neglecting van der Waals effects, the effective 
colloid-colloid interaction takes a hard core Yukawa form 
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where /3 = l/ksT^ Aij = except in the case of non- 
additivity (see below) , Zj^ is the charge on colloid species 
/c, Xb is the Bjerrum length, a^j = (a^ + crj)/2 is the 
mean of the diameters of colloid species z, j and r is the 
centre separation of the two species. The inverse Debye 
screening length is denoted by z^: = ^/^TrXBpion- Linear 
Poisson-Boltzmann theory is only valid for small electro- 
static potentials (colloid charges). However, at higher 
potentials, the pair interaction retains a Yukawa form 
but with a smaller renormalized (effective) charge [7j. 
Using the effective charge, a phase diagram can be calcu- 
lated using Yukawa theory, which accurately reproduces 
experiment, for parameters comparable to ours [8, 9 . 

Treating charged colloids as a one-component system 
with Yukawa repulsions (with fixed parameters), one 
does not expect demixing into colloid-rich and colloid- 
poor phases, as this would typically require some form of 
attraction between the colloids. Experimental observa- 
tions of colloidal fluid-fluid phase separation [10] there- 
fore attracted much attention. Some explanations in- 
voked higher-order terms where the electrostatic repul- 
sions between two colloids might be screened by an inter- 
vening third particle [11 . This 'macroion screening' was 
subsequently experimentally verified [12 , but is a weak 
effect. Other possibilities focussed on one-body terms, 
such as ion condensation, which turn out not to be signif- 
icant for monovalent ions [13] and 'volume' terms, where 
the entropy of small ions, colloid-ion and colloid self in- 
teractions contribute to the free energy [14 . 

Here we consider a mixture of similar sized colloids 
with different charges, which raises the additional pos- 
sibility of non-additivity. Non additivity refers to the 
case where the cross interaction differs from the mean 
interaction between each species, for example in Eq. (IT]) 
Aij 7^ when i j^ j and otherwise. For the nearly size- 
symmetric, dilute colloidal systems we consider, these 
non-additivity effects are expected to be small, so Eq. 
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Figure 1: (color online) (a) Phase boundary in the {(t)H^ps) 
plane for the Hi and H2 one-component systems. Experi- 
mental results for Hi and H2 are denoted by filled circles and 
squares respectively. Predictions from Yukawa theory [18] 
are shown as the solid (Hi) and dashed (H2) lines, (b) Over- 
iew of phase separation {(t)tot = 0.0234, t = 15 hrs). From 
left: X — 0.05 (crystal), 0.30, 0.50 (phase separated) and 
0.95 (fluid), (c) Phase diagram of binary charged colloids de- 
fined by (t)tot and X obtained at t = 240 hrs. Filled squares, 
H-rich crystal (C); open circles, phase separation (demixed); 
filled circles, fluid (F). Dashed line denotes fluid-crystal phase 
boundary in the L-dominated regime. 



^ holds with |Ai^| < 0.01 fUlIIS], and thus we expect 
nearly additive Yukawa behavior. Like one-component 
Yukawa systems, additive binary Yukawa systems are not 
expected to phase separate [TT. Remarkably, our exper- 
iments reveal demixing into phases rich in high-charge 
and low-charge particles respectively. The phase rich in 
high-charge particles is found to be a colloidal crystal. 
By including one-body Wi contributions through volume 
terms O [l^ we identify a driving force for phase sepa- 
ration. 



II. EXPERIMENTAL 

We used low charge colloidal silica (Z = — 170, cr^ = 
110 nm; species L), (Japan Catalyst, Co., Ltd). We syn- 
thesized higher charge polystyrene particles. Hi {Z = 
-870, o-j:/! = 92 nm) and H2 {Z = -940,cr//2 = 90 nm). 
Effective charges were determined from electrical con- 
ductivity measurements as described in [9 . All samples 
were purified by dialysis and ion exchange [19]. Aque- 
ous dispersions of L exhibit no crystallization for volume 



fractions ^^ < 0.07. The phase diagrams for the one- 
component systems of Hi and H2, in the volume fraction 
(I)Hj salt number density ps plane are shown in Fig. CTa) 
[18]. For our particles, any deviation from Eq. CT is 
expected in the high-charge case. Since our results agree 
well with the predictions from Yukawa theory, we con- 
clude that it is reasonable here to treat colloid-colloid in- 
teractions as a Yukawa system with an effective charge. 
No significant difference was seen between Hi and H2 
which are hereafter termed 'H'. 



III. RESULTS 

Phase separation of binary like-charged colloids is 
shown in Fig. [lib). Here the total volume fraction 
^tot = ^L + 4>H = 0.0234. We scan from H-rich to L-rich 
by changing the mixing ratio X = (j)L/(l>tot- When low 
charge particles dominate {X = 0.95) (right), the system 
appears opaque, indicating a colloidal fluid. Meanwhile 
for a majority of high charge particles {X = 0.05) (left), 
the system is iridescent — a colloidal crystal. However at 
X = 0.3, 0.5 we find a macroscopic phase separation into 
iridescent (upper) and disordered (lower) regimes with a 
clear boundary. The change in color indicates that as 
X increases, so does the lattice constant of the colloidal 
crystal. The phase diagram at t = 240 hrs is shown in 
Fig. [ijc). This has a pronounced asymmetry, such that a 
higher (^tot was required for demixing in the case of L-rich 
samples. The phase diagram in Fig. flic) was unchanged 
up to t = 1440 hours. 

We now turn our attention to the composition of each 
regime. Concentration profiles were determined as fol- 
lows for L^H {(l)tot = 0.0234, X = 0.50). By finding 
the concentration of silica monomers, (following decom- 
position with sodium hydroxide) using the molybdenum 
blue method we obtained (J)l{z) (silica) at 15 heights z. 
<I^h{z) (polystyrene) was determined via turbidity mea- 
surements in 90% ethylene glycol whose refractive index 
is matched to silica. Fig. |2la) shows the vertical dis- 
tribution of L (filled circles) and H (open squares) thus 
obtained for the samples at t = 1.5, 60 and 262 hrs. The 
distribution of L becomes inhomogeneous with time due 
to phase separation, while the H particles show some- 
what weaker segregation. At t = 262 hrs, the averaged 
(j)L = 0.05 and (J)h = 0.014 in the upper region while in 
the lower region (J)l = 0.019 and (pn = 0.010. 

We then made a closer investigation of the composi- 
tions of the top and bottom regions, using confocal mi- 
croscopy (Zeiss LSM 510 META). We gelled sucessive 
samples by curing a UV-sensitive additive [20 at times 
t = 1.5, 23 and 60 hrs. The resulting images at heights 
z = 4 and 27 mm reveal the process of phase separa- 
tion for {(l)tot = 0.0234, X = 0.30) [Fig. ^h)]. After 
1.5 hrs, round crystalline domains, a few tens of mi- 
crons in size, are present throughout the sample. In the 
top region of the sample, the average domain size then 
increased to approximately 100 jam up to t = 23 hrs. 
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Figure 2: (color online) (a) Lateral concentration profiles of 
species L and Hi in their binary mixtures {(j)tot — 0.0234, X — 
0.50) at t = 15, 60 and 262 hrs. Filled circles, L; open squares. 
Hi. The profiles for L one-component system are also shown 
for comparison in dashed curves, indicating that sedimen- 
tation is not important on these timescales. (b) Confocal 
images of the system immobilized by polymer gel at differ- 
ent times at two heights. Scale bar applies to all images 
(0^^^ = 0.0234, X = 0.30). Dashed line in (b.v) indicates 
phase boundary. Note that the mixing ratio X is different 
between (a) and (b). (c) Contributions to changes in energy 
density upon phase separation for a mixing ratio X — 0.5. 
Free energy contribution from one-body volume terms 5Wi 
(red circles) dominate over two-body potential energy contri- 
butions 5U2 (pink squares). 



while the lower region became fluid. After 60 hrs, the 
top was mostly crystalline domains. If we suppose that 
the crystals are H-dominated and the fluid which lies be- 
tween is L-dominated, then the local volume fraction of 
highly charged particles in the crystals 0^^ < (J)h{z) the 
mean volume fraction and ^^^ > (t)L{z) where the mean 
(I^l,h{z) are given in Fig. [2|a). This indicates that our 
experiments may be incompletely demixed. The H-rich 
regions have a microstructure of H crystals and L fluid 
in the interstices [Fig|2|b.iii)], while the L-rich fluid also 
contains H particles. 



IV. DISCUSSION 



We now consider a mechanism for the phase separation 
we observe, recalling that such behaviour is not expected 
for an additive binary Yukawa system [17]. Preferential 



sedimentation of the denser silica is ruled out: sedimen- 
tation of a one-component L system occurs on a timescale 
of weeks, as the dashed lines in Fig. [2|a) indicate. By 
contrast phase separation takes hours. We also tuned the 
Z value of the L particles by varying the pH [21^. For 
Zl = 900 we observed no phase separation. We conclude 
that demixing is not related to the different colloid mate- 
rials, but to the difference in charge between the colloids. 

We identify two mechanisms for phase separation: 
crystallisation of one species (a two-body effect) and a 
reduction in potential energy due to ion-colloid coupling 
(a one-body effect). Crystal/crystal and fluid/fluid phase 
separation were not observed. Demixing is then associ- 
ated with crystallisation of the strongly charged species 
H, and we believe this is more than coincidental. Crystal 
nuclei of species H could expel L as the initial phase sep- 
aration mechanism, suggesting some asymmetry in the 
phase diagram as we find [Fig. [ijc)], along with the ab- 
sence of phase separation in the case that neither species 
crystallizes. Charged colloids have strong interaction en- 
ergies and, as mentioned above, additive binary Yukawa 
systems are stable to demixing jlT]. Moreover entropic 
contributions from small ions stabilise the mixed state. 
In other words, without some additional mechanism to 
lower the free energy, we expect that crystallisation of H 
would be arrested by the free energy cost associated with 
phase separation. We therefore enquire as to a mecha- 
nism for reduction in free energy, which would promote 
demixing. 

We proceed by considering the one-body contributions 
to the Hamiltonian Wi . Although these volume terms are 
not related to the colloid coordinates, they capture ionic 
contributions the free energy. In this way, the volume 
terms contribute to equilibria between phases with differ- 
ing compositions. In fact, since the ions are so much more 
numerous than the colloids, they dominate the osmotic 
pressure, and entropic contributions to the free energy, 
while ion-colloid coupling is an important part of the po- 
tential energy [14 . The contribution to the Helmholtz 
free energy density from the volume terms reads [5j 
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where Wi = Vwi and V is the total volume. The salt 
number density ps we take from the self-dissociation of 
water. We compared with higher salt concentrations and 
found very similar results. We use values close to our ex- 
perimental parameters with Zh — —900 and Zi^ = —170, 
and diameters cl = cr^ = 100 nm. The osmotic pres- 
sure is then 11 = '}Zk=i PkPk —Wk where the volume term 
contribution to the chemical potential is /i/^ = dwi/dpk- 
The Grand Canonical formalism behind Eq. ^ in- 
cludes small ion contributions to 11 ^ . Here we assume 
comlete demixing, so equating the osmotic pressures of 
the coexisting phases fixes partitioning of volume, which 
yields the colloid volume fraction of each demixed phase 



0^ and 0^. Equation p| then yields 6wi^ the change 
in free energy density due to the volume terms upon 
phase separation. Swi is shown in Fig. [2|c) for X = 0.5. 
For all colloid volume fractions, Swi is negative, suggest- 
ing that the system has the potential to demix. How- 
ever, the magnitude of Swi increases significantly around 
^tot ^ 0.003, which is similar to the point at which our 
experiments demix [Fig. [Tic)]. 

Now two-body contributions are also density- 
dependent. The prefactor and screening length in 
Eq. (fTl) both depend on the ionic strength, which 
is a function of the colloid volume fraction [22 . We 
therefore also calculate Su2^ the effective colloid-colloid 
contribution to the change in potential energy density, 
using Monte Carlo simulation in the canonical ensemble 
according to Eq. (fTl) [23 . We neglect entropic contri- 
butions to the free energy arising from W2, recalling 
that ionic contributions to Wi are much more signifi- 
cant. The simulations yield the pressure contributions 
from 2-body interactions, but these make almost no 
difference to the values of ^^ and (j)jj at coexistence 
and are neglected. Su2 is plotted in Fig. [2lc). At 



intermediate densities (0.001 



< 



'Hot 



< 



0.01), we find 



evidence that two-body contributions might suppress 
phase separation as Su2 > 0. However at higher volume 
fractions {4>tot ^ 0.01), Su2 is dwarfed by the volume 
term contribution. Our treatment therefore suggests 
a lowering in free energy upon phase separation as a 
similar state point to that found in experiments. 

Considering the terms in Eq. ([2|, we can identify a 
driving force for phase separation. The first term repre- 
sents ion entropy, which tends to oppose phase separa- 
tion. The second term represents the colloid self-energy, 
that is to say ion-colloid interactions. Here, partition- 
ing leads to enhanced ion-colloid coupling in the H-rich 
phase, and we argue that this stablises phase separa- 
tion. There are other, smaller, effects, for example the 
change in Debye screening length in each phase following 
demixing, and its contribution to the potential energy 
contributions to the colloid-colloid terms Su2. We have 
observed that the lattice constant of crystals of H parti- 
cles increases with X in Fig. fllb). This arises naturally 
because increasing X at fixed (j)tot reduces the osmotic 
pressure, allowing the crystals to expand. 

Finally, we consider the generality of our findings. The 
mechanism we propose should be relevant to multicom- 
ponent like-charged systems, such as dusty plasmas and 
protein solutions. Segregation in the latter can have im- 
portant physiological implications. Using parameters of 
aij = 1 nm and net charges of Zh = 6 and Zl = 1 and 
Cs = 0.01 M to describe proteins in physiological condi- 
tions, we find a free energy difference per molecule of 0.1 
ksT at ^tot = 0.13. 



A number of assumptions have been made here, which 
suggest that our findings could be investigated futher in 
the future. In particular, (1) we have assumed complete 
phase separation, (2) that linear Poisson-Boltzmann the- 
ory holds, and (3) that higher order terms such as three- 
body colloid-colloid interactions can be neglected. For 
our parameters, for (ptot ^ 0.01, (3) appears to hold [11 . 
(1) is likely an oversimplification, we expect that some 
mixing of the phases may occur, particularly in the L- 
dominated fluid, however rather than a precise compari- 
son, here we seek to identify a physical driving force for 
phase separation. 



V. CONCLUSION 

In summary, we report that suspensions of colloids of 
differing charge but near-equal size undergo phase sep- 
aration into crystal and fluid phases, which is not ex- 
pected on the basis of the usual Yukawa/DLVO descrip- 
tion. Demixing appears to be initiated by crystallization 
of the strongly charged species and is stablised by one- 
body effects introduced by the small ions. The leading 
driving force is the ion-colloid interaction energy which 
favors partitioning to a higher density of strongly charged 
colloids. 

Apart from our observation that crystallization of one 
species is involved, the mechanism of this novel phase 
separation in binary systems is open to investigation. An 
intriguing question concerns the fact that macroscopic 
phase separation appears to split into partially demixed 
regions of comparable volume: complete phase separa- 
tion should result in a very asymmetric partitioning, to a 
low density phase of high-charge and high density phase 
of low-charge particles. It is possible that this is not 
reached due to kinetic trapping in the experiments. Al- 
though we have offered what we believe to be the most 
compelling explanation, it would be extremely desirable 
to investigate these findings with more sophisticated the- 
oretical treatments. In particular, at the two-body level, 
a full free-energy calculation which includes crystallisa- 
tion of the strongly-charged species would be most desire- 
able. In closing, we hope that the phenomenon of phase 
separation in binary like charged colloids will prove rel- 
evant to understanding of fluid-fluid phase separation in 
charged colloids and to segregation in biological systems 
and dusty plasmas. 
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